Spin and Charge Transport on the Surface of a Topological Insulator 
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We derive diffusion equations, which describe spin-charge coupled transport on the helical metal 
surface of a three-dimensional topological insulator. The main feature of these equations is a large 
magnitude of the spin-charge coupling, which leads to interesting and observable effects. In par- 
ticular, we predict a new magnetoresistance effect, which manifests in a nonohmic correction to a 
voltage drop between a ferromagnetic spin-polarized electrode and a nonmagnetic electrode, placed 
on top of the helical metal. This correction is proportional to the cross-product of the spin polar- 
ization of the ferromagnetic electrode and the charge current between the two electrodes. We also 
demonstrate tunability of this effect by applying a gate voltage, which makes it possible to operate 
the proposed device as a transistor. 



Time-reversal invariant topological insulator (TI) is a 
new state of matter, distinguished from a regular band in- 
sulator by a nontrivial topological invariant, which char- 
acterizes its bandstructure |l[. Its theoretical [2] and 
experimental 3] discovery has accordingly generated a 
great deal of excitement in the condensed matter physics 
community. The most robust observable consequence of 
a nontrivial topological character of these materials is the 
presence of gapless helical edge states, whose gaplessness 
is protected by time-reversal symmetry and is thus ro- 
bust to perturbations that do not break this symmetry. 
In particular, the surface of a three-dimensional (3D) TI, 
such as Bi2Se3 or Bi2Te3 \4\, is a 2D metal, whose band- 
structure consists of an odd number of Dirac cones, cen- 
tered at time-reversal invariant momenta in the surface 
Brillouin zone [2]. Assuming the Fermi surface encloses 
only one Dirac point, the low-energy Hamiltonian, de- 
scribing such a 2D metal, is given by: 
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where z is along the normal to the surface, vp is the Fermi 
velocity, t is the vector of Pauli matrices, summation 
over repeated spin indices is implicit, and we use Ti = 1 
units henceforth. The eigenstates of this Hamiltonian are 
labeled by helicity, i.e. projection of the spin of the elec- 
tron on the direction of its momentum, hence the name 
helical metal (HM). The most obvious physical property 
of a HM is a strong coupling between the spin and orbital 
degrees of freedom, the energy scale characterizing this 
couphng being the Fermi energy ep- Spin-orbit (SO) cou- 
pling of such a magnitude is unprecedented among known 
materials and it is thus extremely interesting to work out 
its possible observable consequences. Some work on this 
subject has already appeared in the literature [a, l6(. 

In this Letter we will focus on spin and charge trans- 
port phenomena in a 2D HM. In this setting, it is use- 
ful to note that the Hamiltonian of a HM Eq.(IT]) is 
very similar to the Rashba Hamiltonian [7|, which has 
been studied extensively in the context of spin trans- 
port in semiconductor-based two-dimensional electron 



gas (2DEG) systems [8l4ld|. The main difference is the 
magnitude of the SO term: while in a semiconductor 
2DEG it is only a weak perturbation on top of the band 
kinetic energy, it is the only term present in the Hamil- 
tonian of the HM. Thus, while the spin-charge coupling 
effects in semiconductor 2DEG's are extremely weak and 
largely experimentally unobservable, they can be ex- 
pected to be significant in a HM. 

We start from the HM Hamiltonian Eq.([T]). We will as- 
sume that the surface Fermi energy, measured relative to 
the Dirac point, is finite (we will take it to be positive for 
concreteness), as is generically the case in real materials 
l] , and that the corresponding 2D Fermi surface encloses 
one Dirac point (we expect our results to hold for any odd 
number of Dirac cones as different cones will contribute 
additively to transport). We will assume that nonmag- 
netic impurities with potential Vi{r) = uo'^a^i^ ~ ^a) 
and 2D density Ui are present at the surface. To derive 
the diffusive transport equation we will adopt the for- 
malism of Ref. 9] , based on the evaluation of the density 
matrix response function. 

The real-time disorder-averaged Green's function of a 
2D HM is a 2 X 2 matrix in spin space. It is convenient to 
write it as a sum of "singlet" and "triplet" contributions 
G^^(k,c^) = Gf-^{k,u:)d^^, + Gf'^(k,c^) • t,,,, where 
R, A stand for retarded and advanced and the singlet and 
triplet Green's functions, projected onto the band with 
positive helicity, are given by: 
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Here ^k = VF\k\ ~ ep and the impurity scattering rate is 
given by 1/r = TmiUQp{ep), where p{ep) — ep/2'irvp 
is the density of states at Fermi energy. We assume 
that epT 3> 1, which allows us to treat disorder per- 
turbatively. We then introduce generalized density op- 
erators gCTCT'(r,i) = *^/(r,i)*g.(r,i), and evaluate the 
retarded density response function X(Tia-2,craai{'r — r',t — 



t') = -i0{t - t'){[g^^,^{r,t),Bl,^{r',t')]). To lead- 
ing order in the small parameter l/epT this is easily 
evaluated by summing all ladder vertex corrections to 
the polarization bubble diagram [9|. Assuming q ^ 
kp and il ^ ep, one obtains the following result 
for the Fourier transformed density response function: 
X{q,n) = x° ~ inTp{eF)I{q,n)V{q,n), where x° is 
the static uniform susceptibility and 1^^10-2, asruiq,^) ~ 

is the diffusion propagator, i.e. the Green's function of 
the spin-charge coupled transport equation we seek to 
derive. Evaluating the momentum integral above, ex- 
panding the result to leading nontrivial order in iflr and 
vp^T, which corresponds physically to coarse-graining 
over length scale of the order of the mean free path 
^ = vpT and time scale of order r, we finally obtain the 
following set of spin-charge coupled transport equations: 
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Here S is the nonequilibrium spin density and N/p(ep) 
is the full local electrochemical potential (we have sub- 
sumed the contribution of the external electrostatic po- 
tential into the definition of the charge density N). D = 
VpT/2 is the diffusion constant and F = ^D/2t = vp/2 
is the spin-charge coupling constant. The out-of-plane 
component of the spin density does not appear in the 
above equations. The reason for this is that any nonequi- 
librium spin polarization in the z-direction will pre- 
cess with frequency of order ep around the momentum- 
dependent SO field z X k and will thus average out to 
zero on time scales longer than 1/ ep. 

It is useful to compare Eq.© with the corresponding 
equations for semiconductor 2DEG systems with Rashba 
SO interactions, derived in Ref.^]. In the latter case, the 
coupling between the in-plane and out-of-plane spin com- 
ponents is large, of order XepT, where A is the Rashba SO 
coupling constant. This leads to interesting and experi- 
mentally observable effects in spin transport [10|. These 
effects will not be observable in HM due to a vanishing 
coupling between in-plane and out-of-plane spin compo- 
nents. In contrast, the spin-charge coupling in Rashba 
2DEG systems is extremely small, with F ^ A(Afci?T)2, 
and is largely experimentally unobservable. In a HM, 
on the other hand, the spin-charge coupling is much 
stronger: V ^ vp, and one can thus expect these effects 
to be easily experimentally accessible. Another differ- 




FIG. 1: (Color online). Schematic picture of the sample. The 
left electrode (FM) is ferromagnetic and can inject a spin- 
polarized current (polarization shown by arrow) into the sam- 
ple. The right electrode (N) is nonmagnetic. The voltage drop 
between the electrodes has a contribution, proportional to the 
cross-product of the polarization and the current. 



ence from the Rashba 2DEG case is an anisotropic na- 
ture of spin diffusion in the HM: the diffusion constant 
is different in the direction of the spin polarization and 
perpendicular to it. This is not unexpected and is also 
a consequence of the strong SO coupling in this system. 
Finally, the relaxation time for the in-plane spin com- 
ponents in the HM case is the same as the momentum 
relaxation time r. In the Rashba 2DEG case, in contrast, 
the relaxation time is much longer, of order t / (Xkpr)"^ . 
Unusual spin-charge transport properties of the HM 
are embodied in a nonstandard relation between the 
charge current and the gradient of the electrochemical po- 
tential, which are no longer simply proportional to each 
other. This can be seen from the first of Eq.Q. Namely, 
since the electron charge is conserved, the charge trans- 
port equation must have the form: 
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where J is the particle current density (equal to 1/ex 
the electrical current density). Then from Eq.(I3]) we can 
immediately read off the following expression for the par- 
ticle current density: 
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The first term here is the usual contribution to the charge 
current, proportional to the charge electrochemical po- 
tential gradient. There is, however, a second term, which 
is proportional to the excess spin electrochemical poten- 
tials, if the spins are polarized in-plane (in the case of a 
clean TI surface this term was first studied in Ref.[5|). 
The physical origin of this term is a definite helicity of 
the Fermi surface states in a HM. 

We will now show that this extra term in the expres- 
sion for the charge current manifests as a novel and po- 
tentially useful magnetoresistance effect. We imagine a 
setup, schematically depicted in Fig.l. Assume that two 
electrodes are deposited on top of a TI sample a distance 
L from each other. Assume the left electrode is ferro- 
magnetic (FM) and can inject a spin-polarized current 



into the surface of the sample, while the right electrode 
is nonmagnetic (N). Let a fixed current / flow between 
the electrodes. We will now use Eq.(I3|) to calculate the 
corresponding voltage drop between the electrodes as a 
function of the degree and direction of the spin polariza- 
tion in the FM electrode. Assuming there is no variation 
in the spin and charge electrochemical potentials in the y- 
direction (i.e. parallel to the electrodes), the spin-charge 
diffusion equations simplify to: 
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We will assume for simplicity that the distance between 
the electrodes L is greater than the spin diffusion length 
in the HM, which is of order £. Then we can assume that 
the spin transmission from the FM to the N electrode is 
zero and take the boundary conditions for Eq.® to be: 
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where J = —DdN/dx — 2TSy is the particle current in 
the x-direction and —{3D/2)dSy/dx is the spin current. 
Note that the term —FN, which, as may naively be con- 
cluded from examining Eq.Q, contributes to the spin 
current, should not in fact be included in the definition 
of the spin current. This term is not zero even in equilib- 
rium and thus represents an equilibrium spin current [ll| , 
which has no relation to measurable transport spin cur- 
rent. One arrives at the same conclusion from a formal 
derivation of the spin current boundary conditions [l^ . 
The second boundary condition assumes that the current, 
injected into the HM from the FM electrode, is spin po- 
larized with the degree of spin polarization parametrized 
by 77. Note that 77 can be either positive or negative, cor- 
responding to spin polarization of the FM electrode along 
y or —y direction. At the N electrode, the spin current is 
taken to be zero, as discussed above. Note that the po- 
larization dependence disappears from Eqs.®, © if the 
FM electrode is polarized in the x-direction, i.e. along 
the direction of the current. 

We now solve Eqs.® with the boundary condi- 
tions (O. The first of Eqs.® simply gives: 
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Expressing dN/dx in terms of 5" and substituting into 
the second Eq.® we obtain: 
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which needs to be solved with boundary conditions ([T]). 
The solution is given by: 
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The first term in Eq.([TO]) corresponds to the spin density, 
injected from the FM electrode. This term decays expo- 
nentially away from the electrode on length scale of the 
order of the mean free path t. The second term, on the 
other hand, is constant and is due to the magnetoelectric 
effect: generation of nonequilibrium in-plane spin polar- 
ization by charge current. Note that the magnetoelectric 
coefficient does not depend on disorder and in this sense 
is universal. 

The voltage drop between the electrodes can now be 
calculated as: 
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Thus the voltage drop consists of two contributions. The 
first contribution corresponds to the usual Ohmic resis- 
tance, which is proportional to the separation between 
the electrodes L and inversely proportional to the con- 
ductivity. The second contribution depends on the degree 
and direction of the polarization of the FM electrode and 
constitutes a new magnetoresistance effect. Note that 
this contribution is proportional to the cross-product of 
the polarization of the FM electrode and the charge cur- 
rent (recall that it disappears when the spin polarization 
of the FM electrode is along the ±x-directions). This 
distinguishes it from all other known magnetoresistance 
effects. The effect has some degree of universality, in 
the sense that it is independent of disorder and the sep- 
aration between the electrodes and only depends on the 
spin polarization of the injected current and the Fermi 
momentum, characterizing the bandstructure of the HM. 
However, the degree of the injected spin polarization of 
course does depend on nonuniversal properties of the in- 
terface between the FM electrode and the HM. 

It is interesting to note that when L < 2£, the resis- 
tance, according to Eq. (fTTI) . may formally become neg- 
ative when 77 is negative. This may mean two things. 
One possibility is that this simply indicates that we are 
going outside of the regime of validity of the gradient ex- 
pansion of the diffusion propagator or of the boundary 
conditions Eq.© at short length scales. Another, more 
interesting possibility is that this signifies a real physical 
phenomenon: a nonequilibrium transition into a zero- 
resistance state, reminiscent of the microwave-induced 
zero resistance states in 2DEG [13| . We leave the resolu- 
tion of this issue to future study. 

It is important to realize that while Ea. ([TT|) was ob- 
tained using a specific form of the spin current boundary 



conditions Eq. ([7]) , our result is in fact largely independent 
of the boundary conditions. If instead of Eq.dT]) we used 
a more general form for the boundary conditions, involv- 
ing both spin densities and spin density gradients (which 
physically would correspond to significant spin-flip scat- 
tering at the electrode- HM interfaces), only numerical co- 
efficient of the magnetoresistance term in Eq. (fTTj) would 
change, while its general form would remain the same. 

The magnetoresistance effect we propose is fundamen- 
tally different from other known magnetoresistance ef- 
fects and is easily distinguishable from them. The extra 
voltage drop, first found by Johnson and Silsbee [14| - 
which always occurs at interfaces between ferromagnetic 
and paramagnetic metals and which results from differ- 
ences in conductivities of the majority and minority spin 
electrons in ferromagnets, is an even function of the spin 
polarization of the ferromagnet (also true when the para- 
magnetic metal is helical), while it is an odd function for 
the effect we propose. The famous giant magnetoresis- 
tance effect [lal requires both electrodes to be FM and 
be separated by a distance, smaller than the spin diffu- 
sion length. The effect we describe here requires only one 
electrode to be FM and does not depend on the distance 
between the electrodes. 

Yet another interesting feature of Eq.(fTT|) is that the 
magnitude of the polarization-dependent term may be 
varied independently of the Ohmic term by applying an 
external gate voltage to the HM to tune the Fermi en- 
ergy. This is based on the observation that the product 
of the Fermi momentum and the mean free path is in 
fact independent of the Fermi energy kpi = 2vp/niUQ, 
which makes the Ohmic term insensitive to changing the 
Fermi energy 10\ . The magnitude of the polarization- 
dependent term, however, does explicitly depend on the 
Fermi energy as '^ l/ep- The possibility of realizing this 
transistor-like effect depends strongly on how insulating 
the bulk TI material actually is. It may thus be diffi- 
cult to realize with currently available 3D TI, in which 
the bulk is not really insulating [l|. We expect, however, 
that these problems will be resolved in the next genera- 
tion of TI materials. 

In conclusion, we have proposed a new magnetoresis- 
tance effect that should occur on the HM surface of a 3D 
TI. The effect manifests itself in a nonohmic correction to 
the voltage drop between two electrodes, placed on top 
of the HM: a polarized FM electrode and a nonmagnetic 
electrode. The correction is proportional to the cross- 
product of the charge current between the two electrodes 
and the spin polarization of the FM electrode and thus 
changes sign when the direction of the spin polarization is 
reversed. We have also proposed a transistor-like device, 
based on this new magnetoresistance effect. 
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